Bosia etal. BMC Systems Biology 201 2, 6:1 31 
http://www.biomedcentral.eom/1 752-0509/6/1 31 



Systems Biology 



RESEARCH ARTICLE OpenAccess 



Gene autoregulation via intronic microRNAs 
and its functions 

Carla Bosia 1 ' 2 * 1 , Matteo Osella 3 ' 4 * 1 , Manama El Baroudi 5 , Davide Cora 2 ' 6 and Michele Caselle 2 ' 7 



Abstract 

Background: MicroRNAs, post-transcriptional repressors of gene expression, play a pivotal role in gene regulatory 
networks. They are involved in core cellular processes and their dysregulation is associated to a broad range of human 
diseases. This paper focus on a minimal microRNA-mediated regulatory circuit, in which a protein-coding gene (host 
gene) is targeted by a microRNA located inside one of its introns. 

Results: Autoregulation via intronic microRNAs is widespread in the human regulatory network, as confirmed by our 
bioinformatic analysis, and can perform several regulatory tasks despite its simple topology. Our analysis, based on 
analytical calculations and simulations, indicates that this circuitry alters the dynamics of the host gene expression, can 
induce complex responses implementing adaptation and Weber's law, and efficiently filters fluctuations propagating 
from the upstream network to the host gene. A fine-tuning of the circuit parameters can optimize each of these 
functions. Interestingly, they are all related to gene expression homeostasis, in agreement with the increasing evidence 
suggesting a role of microRNA regulation in conferring robustness to biological processes. In addition to model 
analysis, we present a list of bioinformatically predicted candidate circuits in human for future experimental tests. 

Conclusions: The results presented here suggest a potentially relevant functional role for negative self-regulation via 
intronic microRNAs, in particular as a homeostatic control mechanism of gene expression. Moreover, the map of 
circuit functions in terms of experimentally measurable parameters, resulting from our analysis, can be a useful 
guideline for possible applications in synthetic biology. 



Background 

microRNAs (miRNAs) are small (about 22 nucleo- 
tides) single-strand RNAs able to interfere post- 
transcriptionally with the protein production of their 
targets. Targeting a vast proportion of protein-coding 
genes [1-3], miRNA-mediated regulation composes an 
important layer in gene regulatory networks. The impli- 
cation of miRNAs in several core cellular processes [4-7] 
as well as in many human diseases [8,9] further confirms 
their biological importance. 

Approximately half of the miRNA genes can be found 
in intergenic regions (between genes), whereas the intra- 
genic miRNAs (inside genes) are predominantly located 
inside introns and usually oriented on the same DNA 
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strand of the host gene [10] (a trend further confirmed by 
our bioinformatic analysis shown in a following section). 
Intergenic miRNA genes present their own promoter 
region [11,12] and their expression is expected to be reg- 
ulated by the same molecular mechanisms that control 
the expression of protein-coding genes. On the other 
hand, experimental and computational results are con- 
sistent with the idea that same-strand intronic miRNAs 
are co-transcribed with their host gene [13-17], and then 
processed to finally become mature functional miRNAs 
[18,19] (although exceptions to this common scheme of 
co-transcripton have been reported [20-22]). 

The host-miRNA co-expression can have a specific 
functional role. In fact, an intronic miRNA can support 
the function of its host gene by silencing genes that are 
functionally antagonistic to the host [23], or more gener- 
ally act synergistically with the host by coordinating the 
expression of genes with related functions [24] . 
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In addition to this "cooperative" miRNA-host rela- 
tion, different studies showed that intronic miRNAs can 
directly regulate the expression of their host gene, estab- 
lishing a negative feedback regulation [10,25,26]. In par- 
ticular, instances of negative autoregulatory feedbacks via 
intronic miRNAs were firstly found by expression analysis 
in human [25]. More recently, two independent large-scale 
bioinformatic analysis, based on different algorithms of 
target prediction, claimed that the occurency of intronic 
miRNA-mediated self-loops (iMSLs) in the human reg- 
ulatory network is significantly higher than expected by 
chance [10,26]. The over-representation of such regula- 
tory module can be interpreted as a sign of evolutionary 
positive selection that has led to an accumulation of a 
specific topology able to perform useful elementary regu- 
latory tasks [27] . In addition, two iMSL circuits have been 
confirmed experimentally: regulation of EGFL7 by its 
intronic miRNA miR-126 [28,29] and regulation of ARPP- 
21 by miR-128b [26]. Both regulations were associated 
to relevant biological functions, the former playing a role 
in cancer proliferation [28], while the latter in vertebrate 
brain physiology [26] . 

The combination of all these pieces of evidence suggests 
that iMSLs are an often exploited and presumably func- 
tionally relevant regulatory circuitry. The open question 
concerns the peculiar functions that an iMSL can accom- 
plish and that could have thus driven their pervasive 
spreading in the human regulatory network. Moreover, it 
would be interesting to understand what specificities of 
post-transcriptional autoregulation by miRNAs can make 
them better suited to fullfil certain tasks with respect to 
the trascriptional self-regulation, so widely used in bacte- 
ria [30]. In this paper we address these questions by mod- 
eling the dynamical and stochastic behaviour of the iMSL 
circuit and comparing its properties to those of alternative 
regulatory strategies such as constitutive expression and 
transcriptional self repression. 

Our results show that, despite of its minimal topol- 
ogy, the iMSL circuitry can implement different biological 
functions. It can speed-up the host gene protein produc- 
tion in response to an activating signal, while delaying its 
switching-off kinetics when the activation drops; it can 
buffer fluctuations coming from the upstream network, 
and generate complex behaviours like a host gene expres- 
sion response obeying " Webers law" (i.e. the magnitude of 
the response depends only on the fold change of the input 
signal). While these different functions can be optimized 
individually, by tuning parameters like molecular pro- 
duction/degradation rates, it will be shown that they all 
represent different ways of making the host gene expres- 
sion robust to external fluctuations. Therefore, autoreg- 
ulation via intronic miRNAs can generally represent an 
efficient homeostatic control of the host gene expres- 
sion, in agreement with the observation that miRNAs are 



often involved in signaling networks to ensure homeosta- 
sis and gene expression robustness [31-34]. In addition 
to model analysis, we present our own bioinformatical 
search for iMSLs in human to further assess their 
statistical over-representation and to propose the best 
predicted candidates for possible future experimental 
tests. 

Besides the understanding of the role of endogenous 
iMSLs, our results can be useful for the growing field of 
synthetic biology [35,36], which has succesfully started 
to make use of RNA-based post-transcriptional regula- 
tions [37,38]. The function- topology map presented in 
this paper can contribute to draw up the manual of bio- 
logical circuits that carry out specific functions for syn- 
thetic engineering, adding a simple and efficient wiring 
strategy that can increase systems' robustness in differ- 
ent conditions. A synthetic realization of an iMSL has 
been indeed recently produced and proven to be effective 
in reducing the expression dependency on gene dosage 
[39]. Therefore, the potential additional functions we will 
show associated to iMSLs could be tested in the near 
future. 

Results and discussion 

Outline of the model 

We are interested in a model of iMSLs that can capture 
the fundamental properties of the circuit, but simpli- 
fied enough to avoid the introduction of too many free 
parameters that would make an exploration of the param- 
eter space unfeasible. In this view, the essential steps of 
transcription, translation, degradation and interactions 
between genes are taken into account as summarized in 
Figure 1A. The host gene is assumed to be under the con- 
trol of an activating transcription factor (TF) with concen- 
tration q, in order to study the dynamical and stochastic 
properties of the circuit in presence of upstream input 
signals. The activation is modeled, as usual in this type 
of descriptions [30,40], representing the transcription rate 
k r (q) of the target as a Michaelis-Menten function of TF 
concentration q: 

k r (q) = (1) 
h r + q 

where k r represents the maximum transcription rate in 
the fully activated state, while h r is a dissociation constant 
specifying the TF concentration at which the transcription 
rate is half of its maximum value. However, the analysis 
can be straightforwardly extended to the case of a Hill 
function (substituting q with q n and h r with h"), if in 
presence of cooperativity. 

On the other hand, miRNAs can exert their action 
repressing translation or inducing degradation of their tar- 
get mRNAs [41]. We construct our model supposing an 
action on target translation. While most of the results 
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Figure 1 Representation of iMSL and the two circuits used for comparison. Schematic views of (A) an intronic miRNA-mediated self-loop 
(iMSL); (B) a gene simply activated by a TF (sTF); (C) a transcriptional self-regulation (tSL). A more detailed representation of the three circuits is on 
the right of the figure. Green rectangles are DNA-genes; s and r are the transcribed miRNAs and mRNAs (orange and red stars respectively) which 
can eventually be degraded (broken grey stars). mRNAs can be translated into proteins p (blue circles) and proteins can be degraded (broken grey 
circles). The reaction rates are reported along the corresponding black arrows: k r (q) and k r (q,p) for transcription; k p (s), k p2 and k p for translation;^, 
g r and g p for degradation. Red arrows represent activations, while red lines ending in bars are repressions. For the iMSL and the sTF, the transcription 
rates are functions of the amount of TFs q, while for the tSL the transcription rate is also a function of the target protein p. In the iMSL, miRNA 
regulation makes the rate of translation a function of the amount of miRNAs s. 



shown in this paper are independent of this choice, some 
dynamical properties of the circuit can actually change if 
miRNA action is mainly due to induction of mRNA degra- 
dation. This issue is discussed in more detail in Additional 
file 1. 

A phenomenological description based on nonlinear 
functions has been proven to be effective in modeling 
RNA interference in mammals [42], and was previously 
applied in computational analysis [43,44]. Along these 



lines, we assume that miRNA regulation makes the target 
translation rate k p (s) a repressive Michaelis-Menten-like 
function of the number of miRNAs (5): 

1 + h 

The dissociation constant h establishes the miRNA level 
at which the target translation rate is half of its maxi- 
mum value k p . With the regulatory interactions denned in 
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Equations 1 and 2, it is possible to represent the dynamics 
of the circuit in Figure 1 A by a set of differential equations: 

dr 

— = k r (q) -g r r 

— = k r (q) -g s s 
dp 

= k p (s) r - g p p, (3) 

where r and p are the levels of host gene mRNA and 
protein products, s is the level of miRNAs, and gi repre- 
sents the degradation rate of the molecular species i. As 
discussed in the introduction, intronic miRNAs (same- 
strand with the host) are expected to be co-transcribed 
with their host gene, therefore their production rate k r (q) 
has the same dependence on the input TF level 

A different representation was introduced in the context 
of bacterial sRNA regulation [45-47] and subsequently 
applied with slight modifications to eukaryotic miRNA 
regulation [48] . In this representation, the degree of catal- 
icity, i.e. the ability of a miRNA to affect multiple mRNAs 
without being degraded, was parametrized explicitly [45] . 
The use of an effective phenomenological function (like 
the one in Equation 2) implicitly assumes a catalytic 
action, as commonly believed for miRNAs [42] . 

Moreover, miRNA biogenesis is a complex and highly 
regulated multi-step process (see [5,18,49] for a review) 
that finally leads to a mature miRNA loaded into the 
RNA-induced silencing complex (RISC) which is actu- 
ally the active form that can downregulate the target 
mRNAs. Recent models of regulation by small RNAs have 
tried to build a more detailed modeling framework taking 
into account these processing steps [39,50], for example 
including miRNA incorporation into the RISC complex 
as a reversible binding reaction. Such detailed model- 
ing approaches have succesfully explained results from 
synthetic biology experiments. However, a more com- 
prehensive mathematical representation of the biological 
processes involved comes at the expense of an increased 
number of free parameters, making more difficult an 
extensive exploration of the parameter space. While in our 
modeling strategy the dynamics of miRNA biogenesis is 
not included, in the supporting information (Additional 
file 1) we report the analysis of the possible effects of the 
delay that the biogenesis process can introduce between 
miRNA transcription and miRNA-mediated repression. 
More generally, the relations between different possible 
models of miRNA regulation are discussed in detail in 
the supporting information (Additional file 1), where it is 
shown that most of the results that will be presented in 
the following are essentially independent on the modeling 
strategy, provided that certain generic conditions on the 
parameters are satisfied. 



In an analogous manner, it is possible to model the 
two circuits that we will use for comparison: a gene 
simply activated by the TF (sTF) without any feed- 
back regulation (scheme in Figure IB) and a transcrip- 
tional self-loop (tSL), in which the negative feedback 
is realized through transcriptional repression (scheme 
in Figure 1C). The properties of each circuit will be 
compared using a so called mathematically controlled 
comparison [30]: all the common parameters will be 
kept to equivalent values, constraining the remainders 
so as to achieve the same steady state of protein 
concentration. 

A deterministic description based on ordinary dif- 
ferential equations can effectively describe the mean 
kinetic behaviour of genetic circuits, thus its predictions 
can be tested with experiments based on averages over 
cell populations. In fact, equivalent mathematical treat- 
ments have correctly predicted the dynamical features 
of several endogenous and synthetics circuits [27,30]. 
However, since gene expression is inherently a stochastic 
process [51-53], we will also make use of a stochastic 
description based on a master equation approach, that 
has Equations 3 as a "mean-field" limit (complete model 
in Additional file 1). To compare the stochastic proper- 
ties and the noise susceptibility of the three regulatory 
strategies in Figure 1, we calculated analytically the rel- 
ative fluctuations in protein level p at steady state and 
confirmed our results with Gillespie simulations (see the 
Methods section for details on simulations). 

Autoregulation via intronic miRNAs has many of 
the structural properties of miRNA-mediated incoherent 
feedforward loops, that represent a diffused and function- 
ally relevant motif in regulatory networks [25,44,54-56]. 
In fact, iMSLs can be thought as a mimimal feedfor- 
wad topology with perfect co-expression of the target 
gene and the miRNA buffering node, and thus, as previ- 
ously observed [39], can be considered a special case of 
post-transcriptional incoherent feedforward loops. As a 
matter of fact, many of the functions that we will show 
the iMSL can perform are consistent with the feedfor- 
ward nature of this circuit, and the analysis presented in 
this paper could be easily generalized to miRNA-mediated 
incoherent feedforward loops, adding new pieces to our 
understanding of microRNA regulation in simple circuits. 

Response times to external signals are altered by 
autoregulation via intronic microRNAs 

The response of a transcriptional unit to a stimulus, such 
as a change in a TF concentration, is steered by the life- 
time of its mRNA and protein products. A fast protein 
turnover speeds up the kinetics, but with a consequent 
high metabolic cost, while in the case of long-living pro- 
teins the timescale of changes in concentration can be 
comparable to the cell cycle time [30,57], which can be of 
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several hours. However, the dynamics of a gene expres- 
sion also depends strongly on the regulatory circuitry in 
which the gene is embedded. For example, it has been 
proven that negative transcriptional self regulation (like 
the one in Figure 1C) and incoherent feedforward loops 
speed up the expression rise-time after induction [57,58], 
while coherent feedforwad loops introduce delays [59]. 

We address in this section the question of how the host 
gene kinetics is changed by being a target of its intronic 
miRNA. To this aim, we consider two opposite simplified 
situations: (i) a sudden activating signal that fully saturates 
the promoter, and (ii) the opposite case of an istantaneous 
drop of the activating signal that completely switches off 
transcription. Case (i) can be studied assuming that at 
t = 0 the transcription rate k r (q) switches from zero to 
its maximum value /c r , and measuring the response time 
Ton defined as the time needed to reach half of the final 
protein steady state. In other words, we integrate numer- 
ically Equations 3 to calculate the time Ton such that 
P(Ton)/Pss = 0.5 (where p ss is the final steady-state pro- 
tein level), starting from the condition r(0) = 5(0) = 
p(0) = 0. In case (ii), in which we assume a drop of 
the activating signal at t = 0, we can similarly define a 
response time Toff looking at the decrease of pit) after 
a switch of the transcription rate from k r to zero at time 
t = 0. The same analysis is performed on a sTF (scheme 
in Figure IB) and a tSL (scheme in Figure 1C) for com- 
parison. The response time To of the simple transcription 
unit sTF is used as a normalization, since Ton(OFF) /To is a 
measure of how much a circuit can alter the response time 
with respect to an unregulated gene. 

Many previous analyses of genetic circuit dynamics have 
assumed short-living mRNAs with respect to proteins 
[57-59]. Within this assumption, the mRNA dynamics can 
be neglected since the timescales are governed by the pro- 
tein kinetics. While this is usually a safe approximation in 
bacteria, in eukaryotes the phenomenology can be more 
complex. In mammals, the mRNA half-life can range from 
minutes to about 24 hours [60,61], with typical values in 
the range of 5 — 10 hours [62,63]. Similarly, protein life- 
times cover quite a wide range, from minutes to several 
days [64]. MiRNAs are usually stable molecules with an 
half-life that can span days [65,66], but there are cases 
of short-living miRNAs, as many miRNAs expressed in 
human brain [67] . Moreover, the miRNA turnover seems 
widely regulated as it happens for mRNAs and proteins 
[68]. In summary, while the situation in which proteins 
are more stable than the corresponding transcripts could 
still be frequent, a variety of specific cases is expected. 
Therefore, we decided to take into account the mRNA 
dynamics and explore different regimes of molecules' half- 
lives. Indeed, we will show that the dynamical response of 
the iMSL circuit depends crucially on the ratio between 
mRNA and miRNA half-lives (r r /rs). In Figure 2A, the 



normalized response time Ton I To to activation is plotted 
as a function of the repression level measured as p ss /po> 
where p ss is the final steady-state protein concentration, 
while po is the steady-state protein concentration in 
absence of negative regulation. The response time of the 
iMSL (continuous lines) and the tSL (dashed lines) is 
reported for different values of the half-life ratio r r /r 5 . As 
a first result, the iMSL can speed up the response time 
with a comparable efficiency with respect to its transcrip- 
tional counterpart, especially when mRNAs are degraded 
fastly enough. On the other hand, when miRNAs are 
short-living with respect to mRNAs, they will reach their 
final concentration faster than mRNAs, thus blocking 
more quickly the initial rise in target protein concen- 
tration. Therefore, the timescales of mRNA and miRNA 
dynamics, determined by their half-lives, define the cir- 
cuit performance in speeding up the response, as reported 
in Figure 2A. In Figure 2C an example of the dynamics 
is reported, showing an acceleration of the response for 
both self-regulation strategies at an intermediate level of 
repression. 

As the repression increases, the response acceleration 
to an activating signal relies more and more on an over- 
shoot of protein concentration, well above the final steady 
state, both for iMSLs and tSLs. If the input signal has to 
drive the host gene to its functional steady state, a large 
overshoot can be unwanted since it represents an unnec- 
essary metabolic cost and a possible source of toxic effects 
[57]. Thus, there is probably a limitation in the repres- 
sion strength that can be applied to minimize the time 
separation between two functional steady states. On the 
other hand, a regime of strong repression makes the iMSL 
a pulse generator, a dynamical feature that can eventu- 
ally lead to adaptation as will be discussed in a following 
section. The observation that the iMSL can function as a 
pulse generator is fully consistent with the fact that this 
circuit is a particular case of incoherent feedforwad loops 
that were indeed previously associated to pulse generation 
[30]. 

While the speeding up of activation is a property shared 
by incoherent feedforward loops (and specifically by 
iMSLs) and tSLs, an interesting peculiarity of the miRNA- 
mediated regulation in iMSLs emerges looking at the time 
required for p concentration to reach zero, starting from 
a constitutive level (Figure 2D reports an example of this 
dynamics). The iMSL can delay the switch-off kinetics of 
the host in the same repression regime where it can accel- 
erate the activation and the extent of the introduced delay 
is again dependent on the mRNA to miRNA lifetime ratio 
(Figure 2B). This apparently counterintuitive behaviour 
can be easily qualitatively explained. When a constitu- 
tively expressed gene senses a transcription stop signal, 
the velocity of protein concentration decrease is estab- 
lished only by protein and mRNA degradation rates. For 
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Figure 2 Autoregulation via intronic miRNAs speeds up the host gene activation and delays its deactivation. (A) Activation response time: 
response time Ton, normalized by sTF response time 7b, plotted as a function of repression level p ss /po {p ss /Po = 1 means no repression) for 
different values of mRNA/miRNA lifetimes. Both iMSL and tSL (continuous and dashed lines) accelerate the response time with respect to sTF (blue 
line). Colors correspond to values of mRNA and miRNA relative stability (r r /r s ), achieved keeping mRNA and protein degradation rates {g r , g p ) fixed 
while varying miRNA degradation rate {g s ). The half-life ratio r r /r s affects also the tSL dynamics, but the dependence is weak (the corresponding 
curves tend to collapse). (B) Deactivation response times: response time Toff, normalized by To, is shown for different repression levels. The blue 
line corresponds to the sTF.The response time for iMSL is plotted for different miRNA half-lives (same color code of A). The iMSL induces a delayed 
host gene response in the same repression range that consents a Ton reduction. (C) Example of target protein temporal evolution in activation for 
the three circuitries. The parameter values correspond to stars in plot A. Time is in protein half-life units. (D) Example of iMSL target protein 
dynamics in deactivation for different mRNA/miRNA half-lives (stars in plot B). The parameter setting for this panel is the following: protein half-life 
t p = 8 hours, mRNA half-life r r = 30 minutes, h = 1 000, k r = 0.21 281 9 s~ 1 , k p = 0.0048 s~ 1 . 



example, long living mRNAs are more persistent and can 
be translated for a longer time after the stop of transcrip- 
tion, and long living proteins are obviously more resilient. 
The same is true for tSLs or transcriptional feedforward 
loops: as the transcription is switched off, transcriptional 
repressors cannot exert any regulation and the protein 
level simply undergoes the exponential decrease dictated 
by mRNA and protein degradation. On the other hand, 
thanks to the post-transcriptional regulation in iMSLs, 
for each single miRNA that is degraded the still present 
mRNAs sense an increase in their translation rate. This 
increase clearly depends on the repression strength that 
miRNAs can exert (thus on the repression fold p ss /po) and 
on the relative stability of mRNAs and miRNAs (r r /r 5 ), as 
a fast miRNA turnover leads to a higher translation rate of 
the remaining mRNAs. Eventually, the general increase in 
mRNA translation rate for each miRNA degradation event 
can lead to a temporary boost in protein concentration 
above the original steady state (see Figure 2D). 



It is important to notice that the dynamics just described 
can be altered if the miRNA acts mostly on mRNA degra- 
dation and depends on the timescale of miRNA-mRNA 
binding-unbinding. While the iMSL can always speed up 
the host gene expression in activation, the delay in the 
switch-off dynamics can vanish in case of fast miRNA- 
mediated induction of mRNA degradation. This issue is 
discussed in more details in the Additional file 1. 

The circuit response dynamics can robustly keep the host 
gene in a high-expression state 

In the regime of comparable mRNA and miRNA life- 
times (red curves in Figure 2) the iMSL circuit can both 
accelerate the response to a switch-on signal and delay 
the switch-off kinetics. This alteration of the dynamics 
makes the host gene ON-state (expression at maximum 
rate) robust with respect to a transient fading of the input 
activating signal, as the one depicted in Figure 3A. In 
fact, the response to an input fluctuation toward zero 
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Figure 3 miRNA-mediated self-loops can keep the host gene 
expression robustly in a ON-state. (A) Schematic representation of 
a transient drop of the input TF q of duration T* (time is in protein 
half-life units r p ) (B) Response of the three circuits to the temporary 
absence of signal depicted in A. The iMSL response (red line) is a slow 
protein concentration decrease, followed by a quick recovery of the 
ON-steady-state when the fluctuation is over. For a tSL (orange line) 
or a sTF (blue line) the switch-off dynamics is just due to mRNA and 
protein degradation. Even if the transcriptional negative feedback can 
accelerate the recovery, the distance d from the ON steady state that 
is reached by the target protein level during the temporary absence 
of the input activator is determined by the switch-off response. 
(C) The distance d from the ON steady-state reached by the target 
level is plotted as a function of the duration of theTF q absence 
(d = 1 when the protein level reaches zero). The iMSL circuit (red line) 
requires a more persistent absence of signal to show a significative 
reduction of the host protein product level with respect to the tSL or 
the sTF (blue line). The parameter values are the same of Figure 2, 
with comparable mRNA and miRNA stability (r r /r s = 1). 

is a slow protein concentration decrease, followed by a 
quick recovery of the ON-steady-state when the fluctua- 
tion is over (Figure 3B). Only a persistent absence of signal 
would cause a complete disappearance of the host protein 
product. In this way, the cell could prevent a drop in con- 
centration of physiologically necessary proteins in merely 
presence of activator fluctuations. A resilient ON-state 
can be biologically important if it ensures the homeo- 
static protein level that must be robustly kept to mantain 
the correct phenotype or if the deactivation/reactivation 



is a costy process that have to be engaged only when 
undoubtedly necessary. 

This property can be measured more quantitatively by 
the distance d from the ON-steady-state that is reached by 
the target protein level during a temporary absence of the 
input activator lasting a time T*. As shown in Figure 3C, 
the iMSL regulation keeps the host gene protein product 
close to its steady state in presence of input fluctuations 
that would almost switch-off a gene transcriptionally self- 
regulated or constitutively expressed. 

Intronic microRNAs, targeting their host gene, can 
implement adaptation and Weber's law 
Adaptation 

Adaptation is defined as the ability of a system to respond 
to a change in the input but subsequently return to the 
original level, even if the stimulus persists. Adaptation is 
ubiquitous in signaling systems. Examples of nearly per- 
fect adaptation range from chemotaxis in bacteria [69] to 
sensor cells in higher organisms [70]. In all these systems, 
the benefit of adaptation can be summarized as the pos- 
sibility of signal detection irrespective of the background 
level, thus widening the range of accessible signals and 
keeping the system robust in presence of perturbations. 

Simple network topologies, as negative feedback loops 
with a buffering node or incoherent feedforward loops, 
can be at the basis of the cellular implementation of adap- 
tation [71]. As a special case of incoherent feedforward 
loops, also iMSLs are expected to be suitable to imple- 
ment adaptation. In this section, we investigate whether 
and in what conditions the minimal topology of a post- 
transcriptional self-regulation through intronic miRNAs 
can perform adaptation. 

It is easy to show analytically (see Additional file 1) 
that in the regime of strong repression (s/h ^> 1 in the 
Michaelis-Menten function in Equation 2) the steady state 
of p concentration is independent of the input level q, 
which is clearly a hallmark of perfect adaptation [72]: after 
an eventual dynamical response to a change in q, the sys- 
tem always returns to its original equilibrium level. On 
the other hand, it is impossible to achieve such an inde- 
pendence on the input level at equilibrium using a tSL, as 
confirmed by the fact that in general circuits with just two 
molecular species are not adaptive [71]. 

More generally, we can evaluate the efficiency in per- 
forming adaptation giving the circuit a step function as 
input and calculating the two indexes of precision P and 
sensitivity S [71,73] represented in Figure 4A and defined 
by: 

p = (Pi -Po)/Po ~ l 
(qi - qo)/qo 

g _ Pmax ~ PO ^ 
PO 
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Figure 4 Adaptation and Weber's law implementation via intronic miRNA-mediated self-regulation. (A) Schematic view of adaptative 
behaviour for a step-like input. S and P are the sensitivity and precision measures described in the main text. (B) Schematic view of Weber's law 
implementation for a two-step input function. E is the error in fold-change detection, as defined in the main text. (C) A summarizing heat-map of 
the iMSL performances in implementing adaptation and Weber's law, as a function of the effective activation q/h r and the repression strength 1 /h. 
The grey region is the adaptive region (P > 1 0 and 5 > 2CV P ), while the region where the system implements also Weber's law (£ < 0.1 ) is depicted 
with a color code representing the E value as reported in the legend. In the red zone the system is not sensitive enough to input variations 
(5 < 2CV P ). On the right, the target protein-level response to a two-step input is reported for the parameters values identified by the corresponding 
lower-case letters in the heat-map. The shaded regions correspond to the 2CV P sensitivity threshold, showing that for a too strong repression the 
circuit response cannot produce a signal beyond the noise level (plots in the red rectangle). The parameter setting is the following: mRNA and 
protein half-lifes as in Figure 2, miRNA half-life is r s = 8 hours, k r = 2.1 281 9 s _1 , k p = 0.048 s _1 , the input function starts from an initial value of 
q 0 = 40 and makes two consecutive steps with fold-change F = 4. 



P is a measure of the difference in the steady-state lev- 
els before and after the stimulus, therefore it is actually an 
estimate of the degree of adaptation. Following [71], we 
define the minimal threshold P > 10 to select adaptive cir- 
cuits. A high value of P is not enough to define adaptation 
since it could merely be a consequence of complete insen- 
sitivity to input changes. Thus, it is necessary to check if 
the peak in pit) concentration is an effective recognizable 
signal. This condition can be formalized requiring a sensi- 
tivity S above the noise level {CV p = a p /(/?)) oip at steady 
state, as can be calculated using the stochastic version of 
the model (see Additional file 1). In particular, we choose 
the threshold S > 2CV p (assuming a noise in the input 
level CV q = 10%) to define a circuit "sensitive" to the input 
signal. 

The general requirement of strong repression for an 
effective implementation of an adaptive response is in 



agreement with the results of an analogous theoreti- 
cal analysis previously performed with detailed models 
of different feedforward loop circuits, including iMSLs 
[39]. In particular, the prediction that stronger repression 
should lead to better adaptivity was tested experimentally 
with a synthetic transcriptional feedforward loop, using 
the number of DNA molecules coding for the circuit to 
modulate the input signal [39] . 

Weber's law 

Certain adaptive systems, besides the ability to return to 
their original value after a signal response, present also 
a degree of response that is proportional to the relative 
change in the input signal and not to its absolute value. 
This feature is known as Weber's law, originally intro- 
duced in the context of human sensory response. Recently, 
this dependence on input fold-change was demonstrated 
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experimentally in eukaryotic signaling systems [74,75], 
and theoretically the feedforward loop topology was 
proposed as a candidate to Webers law implementation in 
gene regulatory networks [76] . Once again, it is natural to 
examine in what conditions also the minimal iMSL circuit 
can satisfy Webers law. 

It is possible to show analytically (see Additional file 1) 
that iMSLs are responsive to input fold-change if three 
conditions are satisfied: 

• Strong repression: s/h ^> 1 k p (s) ~ k p h/s 
(condition for perfect adaptation), 

• Almost linear promoter activation k r (q) ^ qk r /h r , 

• Fast mRNA dynamics (short mRNA half-life with 
respect to miRNA and protein ones): r(t) —> r ss . 

As for the case of adaptation, we can quantify the effi- 
ciency in Webers law implementation for a generic set 
of biochemical parameters. To this aim, a two step input 
function is provided such that each step has the same fold- 
change but different background levels (see Figure 4B). 
As previously proposed [76], the error E in recognition of 
fold changes can be quantified using the difference in the 
response peaks: 

j-, Pmax2 ~ Pmaxl /r x 
t = . (b) 

Pmaxl 

Parameter space of adaptation and Weber's law 

Using the observables denned in Equations 4 and 5, it 
is possible to explore the conditions in which adaptation 
and Webers law are successfully performed by iMSLs. 
An illustrative example is depicted in Figure 4C, where 
two effective parameters are varied: the effective promoter 
activation q/h r , and 1/h which measures the repression 
strength since h is the number of miRNAs necessary to 
reduce to one half the target translation rate. The grey 
region depicts the parameter space where precise adap- 
tation is performed (P > 10), while in the excluded red 
region the dynamical response of the circuit is not able 
to go beyond the noise level (S < 2CV p ). The E value is 
reported with the color code in the legend when the min- 
imal condition E < 0.1 is satisfied, i.e the two steps of the 
input produce the same response within 10%. Adaptation 
and Weber's law can be encoded by iMSLs in a param- 
eter region that span several orders of magnitude of the 
effective parameters. Therefore, the only constraint is that 
the effective parameters have to approach the appropriate 
limits, without the need of fine-tuning. 

It should be noticed that the general condition of strong 
repression required for both functions is limited by the 
circuit sensitivity. This is partially due to the fact that a 
too strong repression can rise the noise level of the cir- 
cuit (see next section) making the achievement of a signal 
significantly above fluctuations harder. 



It is interesting to consider what are the functional 
advantages that these two functions can provide to the 
host gene. Both adaptation and Webers law can bestow 
robustness to the expression program of the host gene. 
An expression state that is not influenced by constant 
inputs thanks to adaptation is robust with respect to the 
ubiquitous cell-to-cell variability in TF concentrations, 
but it is still responsive to signals that induce dynam- 
ical variations of TF levels. When additionally Webers 
law is implemented, also the dynamical response can be 
kept homogeneous in a cell population. In fact, in this 
case the response profile is only due to the input fold- 
change and not on its absolute value that is affected by 
the potentially variable background level [76] . Moreover, 
Webers law naturally encodes a noise filter. In fact, since 
the noise level is expected to scale with the background TF 
concentration, a dependence on fold-change can rescale 
appropriately the threshold at which the response is trig- 
gered, thus allowing a better signal/noise discrimination 
in different background conditions [76]. 

Autoregulation via intronic microRNAs reduces the host 
gene expression fluctuations 

All the functions of iMSLs discussed so far contribute to 
enhance the robustness of the host gene expression. It is 
therefore natural to analyze a stochastic model of iMSLs 
to test directly their ability to filter out fluctuations. The 
stochastic analysis of the system is reported in detail in 
Additional file 1. The results in terms of noise-buffering 
properties at the steady state for the iMSL are similar to 
those obtained for the incoherent miRNA-mediated feed- 
forward loops (see [44]). By filtering fluctuations prop- 
agating from the upstream TF, the steady-state target 
protein level achieved with an iMSL is less noisy than the 
same target amount obtained with a simple sTF or a tSL 
(Figure 5A,B). In particular, the target noise CV p for the 
iMSL shows a U-shaped profile with a well defined min- 
imum, thus allowing us to identify the parameter values 
that optimize the noise reduction properties (Figure 5B). 
This prediction could in principle be tested tuning the 
repression strength, as shown in [77] for a tSL. Also a 
tSL can in fact optimally filter noise for well defined val- 
ues of repression strength [77-79], as shown in Figure 5B 
(orange dots and line). For this circuit the mechanism is 
well understood: an excessive increase of the repression 
strength (while potentially improving the noise reduction 
of the circuit) reduces the copy number of mRNAs and 
proteins with a consequent rise in intrinsic fluctuations 
(which can overcome attenuation). Thus, there is just a 
well denned range of repression strength for which the 
noise reduction is optimal, as shown in experiments [77]. 

It is interesting to notice that, even if iMSLs and tSLs 
show similar noise reduction properties, the miRNA- 
mediated self-regulation actually performs better than 
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Figure 5 Intronic miRNAs can buffer noise in host gene expression. (A) An example of the target protein distributions for the three circuits 
(repression level (p)/(po) = 0.2). Lines are gamma distributions with first two moments calculated analytically, while histograms are the result of 
Gillespie simulations. The distribution for the iMSL circuit (red line and histogram) is the narrowest, showing that, even if also a tSL (orange line and 
histogram) can reduce noise with respect to a sTF (blue line and histogram), the iMSL is outperforming. (B) Target noise CV P as a function of the 
repression level (p)/(po) for the three circuits. Lines are analytical predictions, while dots are the result of Gillespie simulations. Given a noise level 
CV q ~ 0.2 in the upstream transcription factor, both the iMSL (red lines and dots) and the tSL (orange line and dots) shows a minimum of noise 
reduction with respect to the sTF (blue line and dots), but the level of fluctuations in the iMSL case is clearly lower. (C,D) Noise reduction on the 
target protein level achieved by the iMSL and the tSL respectively. The noise reduction CV p /CV s p TF (where CV s p TF measures the fluctuations around 
the same mean level for a sTF) is evaluated at different degrees of transcriptional activation {q)/h r and repression (p)/(p 0 ).The same color gradient 
is used in both heat maps, showing that the iMSL reduces fluctuations on a larger parameter region and to a greater extent. In the white regions 
CV P > CV S p F .JUe parameter values are the following: mRNA half-life r r = r w = 30 minutes, protein half-life r p = r q = 1.5 hours, k w = 3.4 10 _3 s _1 , 
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(see Additional file 1 for the definition of kw and kq). 



the transcriptional self-regulation. As it is possible to 
see in Figure 5A (where histograms and continuous 
lines are respectively the result of Gillespie simulations 
with full nonlinear dynamics and gamma distributions 
with analytically calculated moments), the probabil- 
ity distributions of the target protein level for the 
three circuits are different. Both autoregulatory cir- 
cuits lead to a target distribution less sparse than 
a sTF, showing that they effectively reduce fluctua- 
tions, but the iMSL distribution is clearly more peaked 
than the tSL one. Similarly, both self-regulation strate- 
gies show an optimal noise buffering for an interme- 
diate repression strength, but again the attenuation is 
larger in the miRNA-mediated case (see Figure 5B). 



This is more clearly shown in Figure 5C,D, where the 
noise reduction CV p /CV^- F (with CV^ TF representing 
the target noise in the case of a simple transcription unit 
producing the same mean amount of proteins) is reported 
for the two autoregulatory circuits. Noise reduction is 
explored for different levels of transcriptional activation 
({q)/h r ) and target repression ((p)/{po), where (po) is 
the target mean value in absence of repression) to shed 
light into noise control and target suppression interde- 
pendence. Averages are here intended at steady state, thus 
the repression level measure (p)/(po) for the stochas- 
tic model is perfectly equivalent to the one used in the 
response time analysis (see Figure 2). In the regime where 
the target is more sensitive to TF fluctuations, i.e. q is 
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far from saturating the promoter, the iMSL can reduce 
the fluctuations up to a factor 0.5 (Figure 5C), while the 
tSL (Figure 5D) is much less effective. Moreover, the heat 
maps in Figures 5C,D indicate that the iMSL can buffer 
fluctuations over a wider range of conditions as well as to 
a greater extent. 

As pointed out in [44], an optimal miRNA-mediated 
noise buffering does not necessarily require a strong 
repression. Indeed, Figure 5C shows that a reduction of 
the mean protein expression to 50% of its constitutive 
level is sufficient to reduce the noise by approximately 
40%. This means that the intronic miRNA can keep the 
expression of its host gene in its homeostatic regime, 
while filtering out fluctuations, without exerting a strong 
reduction of its concentration. This result agrees well with 
the observation that miRNAs act often to fine-tune their 
targets rather than to switch them off completely [80] . 

Sketch of the one-to-many topology-function map 

This section summarizes the functions found to be asso- 
ciated to intronic miRNA-mediated self-loops into a qual- 
itative "map of functions", showing the different, although 
overlapping, ranges of biochemical parameters in which 
each specific function is optimized. The emerging map 
between parameter values and functions can be useful to 
understand the presence of the iMSL architecture in dif- 
ferent biological contexts and gives general guidelines for 
the design of synthetic circuits with a desired behaviour, 
well beyond the simple suggestion of a network topology. 

As Figure 6A shows, strong repression ((p)/(po) 1) 
is a general requirement for the implementation of adap- 
tation and Webers law, the latter additionally requiring 
an almost linear activation of transcription ((q)/h r 1). 
A sufficiently strong repression is also required to confer 
robustness to the high-expression state (induced by strong 
activation {q)/h r ^> 1) of the host gene in presence of 
input temporary drops. On the other hand, for intermedi- 
ate host activation, where the host gene promoter is highly 
sensitive to changes in the TF concentration, the iMSL can 
efficiently buffer fluctuations at steady state without the 
need of strong repression. 

Looking at a finer scale the strong repression regime 
((p)/(po) < 1/2), a smooth transition in the dynamical 
behaviour of the circuit can be observed (see Figure 6B). 
At first, the host gene is able to fastly transit between 
two well distinct steady states after induction. When the 
repression is further increased, this fast ON-activation 
relies increasingly on an overshoot well above the final 
equilibrium at which the dynamics asympotically relaxes. 
Therefore, the concentration profile resembles a pulse. 
Finally, for high enough repression the system returns 
to the initial steady state after the pulse, a necessary 
condition for the implementation of adaptation and 
Webers law. 



The relative half-life of the molecules involved, in par- 
ticular of miRNAs and mRNAs, is another ingredient 
that can strongly influence the dynamical behaviour (see 
Figure 6C). For example, a miRNA half-life comparable 
to the mRNA one allows a trade-off between accelera- 
tion of the ON-dynamics and delay of the OFF-dynamics, 
making the state of high expression of the host gene 
robust to fluctuations. On the other hand, mRNA life- 
time must be short with respect to the other molecules 
lifetimes for a dynamical response following Webers 
law. 

The present analysis of the iMSL functions considers the 
circuit as isolated, while realistically a single microRNA 
can target hundreds of genes. As recently pointed out, 
the degree of repression of a target depends on the level 
of expression of all possible target genes [45,81], since 
their mRNAs can dilute the pool of available miRNAs. 
Therefore, the expression profile of alternative miRNA 
targets is a variable that can potentially alter the dynam- 
ics of iMSLs (as shown for incoherent feedforward loops 
[44]), and thus have to be carefully taken into account in 
experimental tests on endogenous iMSLs. 

Identification of intronic miRNA-mediated self-loops in 
human 

In this section, we briefly describe our bioinformatic 
search of iMSLs. Our main goal is to provide an updated 
list of candidates to eventually test our theoretical pre- 
dictions. We performed a genome wide search of intronic 
miRNA-mediated self-loops along the lines of two papers 
which recently addressed the same issue [10,26]. The dif- 
ferences between our results and those quoted in [10,26] 
are mainly due to a different choice of the algorithms used 
to predict miRNA targets and in some cases to the use 
of updated versions of the corresponding databases. We 
identified the same strand intronic miRNAs using as ref- 
erence the Ensembl (release 57) database. A summary of 
our results is reported in Table SI of Additional file 1 and 
in Figure 7A, where the percentage of intergenic versus 
intragenic miRNAs is plotted and, for the intragenic ones, 
the relative ratio of exonic versus intronic miRNAs and 
of same-strand versus opposite strand is also reported. 
Target identification was performed using 8 different algo- 
rithms: TargetScan human v. 5.0 [3,82], miRanda - release 
2008 [83,84], RNA22 [85], PITA-4way [86], MirTarget2 
[87], PicTar [88], Diana microT v.3 [89] and TargetMiner 
v.l [90]. In this respects, our analysis could be considered 
as a combination and extension of the one reported in ref- 
erence [10] (where only the first 6 algorithms were used) 
and the study of [26] (where the authors used their own 
prediction algorithm). We were able in this way to find 
a total of 77 iMSLs confirmed by at least one algorithm 
(details are reported in Table S2 of Additional file 1). Since 
these algorithms are very different, we did not try to give 
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Figure 6 Map of functions for an intronic miRNA-mediated self-loop. (A) An ON-state of host gene expression is defined by full promoter 
induction {{q)/h r » 1), and a sufficiently strong miRNA repression {(p)/(po) < 0.5) can keep it robust in presence of temporary drops in the 
activator concentration. In the strong repression regime ((p)/(po) < 1 ) adaptation can be observed, and for almost linear transcriptional activation 
{{q)/h r < 1) the host gene response can show an adaptive dynamics following Weber's law. Fluctuations can propagate from the upstream TF 
more efficiently if the target promoter is highly sensitive to changes in TF level {(q)/h r % 1), thus in this parameter region noise buffering is more 
relevant, with a maximum in efficiency for intermediate repression ((p)/(po) ~ 0-3). (B) A zoom on the strong repression region shows a transition 
between different dynamics. A step input can induce a fast transition of the host gene expression between two distinct steady states, but increasing 
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an absolute score to our results, but ordered them start- 
ing from those which were assessed by the largest number 
of target prediction methods (Figure 7B and Table S2 of 
Additional file 1). Following a standard recipe (see [10] 
for a similar choice) we consider the number of differ- 
ent algorithms that agree on a certain target prediction as 
a measure of the confidence of such prediction. Interest- 
ingly, 28 of our iMSL agree with previous predictions of 
iMSLs [26] and for two of them an experimental valida- 
tion of the miRNA-host gene regulation exist [26,28,29]. 
Moreover, a recent study [17] provides evidence support- 
ing a feedback mechanism between miR-438 and IGF2, 



in agreement with our list of iMSLs predicted by only 
one method. In order to test if these iMSLs are over- 
represented we performed two independent enrichment 
tests. First we performed a reshuffling of the host genes 
while keeping miRNA target predictions unchanged. Sec- 
ond, we randomized the union of the datasets of miRNA 
target predictions obtained with the eight algorithms dis- 
cussed above, keeping the host genes unchanged. In both 
cases we evaluated the Z score which turned out to be 
Z = 4.63 for the first test and Z = 5.52 for the sec- 
ond one (see the Methods section for more details). The 
results of the tests are plotted in Figure 7C. They suggest, 
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Figure 7 Classification of miRNAs and randomization results. (A) Classification of human miRNAs based on the percentage of intergenic and 
intragenic miRNAs. Intragenic miRNAs are divided in exonic and intronic miRNAs and each group can be further classified as same strand or opposite 
strand. All UTR miRNAs were included in the group of the exonic miRNAs. Data are reported in Table SI of Additional file 1 . (B) Number of intronic 
miRNA-mediated self-loops as a function of the number of target prediction methods in agreement: we found a total of 77 iMSLs predicted by at least 
one prediction methods, 25 of them are predicted by at least two different methods and only three of them are predicted by 5 different methods. 
(C) Results of the permutation test: the number of iMSLs in the human network is plotted as a dashed line alongside the distributions (normalized 
histograms) of the number of iMSLs found using the two randomization strategies (described in the main text) over 1000 experiment repetitions. 



in agreement with what already observed in [20,26], that 
this particular class of network motifs is under positive 
selection. 

Conclusions 

This study presents a fairly comprehensive survey of the 
possibile functions associated to a miRNA-mediated cir- 
cuit composed of one protein-coding gene (host gene) 
negatively regulated by a miRNA located in one of its 
introns. In particular, we have shown that, thanks to the 
miRNA-mediated self-regulation, the host gene expres- 
sion responds to input changes with an altered timing, 
its response can be adaptative and follow Weber's law, 
and fluctuations propagating from the upstream network 
are buffered at steady state. Each of these functions can 
confer robustness to the expression program of the host 
gene, suggesting that miRNA-mediated self-loops repre- 
sents a simple homeostatic control. For example, adapta- 
tion makes the host gene expression level at equilibrium 
independent of the cell-to-cell variability of transcription 
factor expression without compromising its sensitivity to 
input changes. Similarly, the host expression dynamics, as 
modified by miRNA autorepression, can mantain the host 



gene in a high-expression state in the face of downwards 
fluctuations in activators' concentration. The association 
of miRNA-mediated self-loops with functions with dif- 
ferent specificities, but apparently the same final aim, 
suggests that, depending on the desired level of the host 
gene expression and on the type of fluctuations that have 
to be more frequently filtered out, the details of the regula- 
tory interactions and the characteristics of the molecules 
involved could have been fine-tuned over evolutionary 
timescales accordingly. Such a fine-tuning of expression 
parameters has been shown to be possible even over short 
timescales in in vitro evolutionary experiments [91]. 

The comparison with an unregulated transcriptional 
unit and with a transcriptional negative feedback indi- 
cates that the specificities of miRNA regulation makes the 
post-transcriptional circuit better suited to implement a 
homeostatic control. This result is in line with the accu- 
mulating clues that miRNAs can help the cell to function 
reliably in presence of perturbations [31-34]. 

Finally, our systematic analysis of the constraints on 
biochemical parameters necessary to optimize each func- 
tion can guide the realization of synthetic versions of 
miRNA-mediated self-loops, as well as contribute to the 
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understanding of the role of their many occurrences in 
endogenous networks. In this perspective, we also pro- 
vide a list of bioinformatically predicted miRN A- mediated 
self-loops in human for future experimental tests. 

Methods 

Stochastic simulations 

Simulations were implemented by using Gillespie's first 
reaction algorithm [92], The reactions simulated are 
those presented in Figure 1 with additional transcrip- 
tion, translation and degradations for the input transcrip- 
tion factor q. Reactions that depend on a regulator were 
allowed to have as rates the corresponding full nonlin- 
ear functions. Results in Figure 5 are at steady state, 
which is assumed to be reached when the determinis- 
tic evolution of the system in analysis is at a distance 
from the steady state (its asymptotic value) smaller than 
its 0.05% (more than 10 times the protein half-life). 
Each data point or histogram is the result of 100000 
trials. 



Bioinformatic methods 

In order to identify human intragenic miRNAs, and asso- 
ciate them to their host genes, we used Ensembl-release 
57 database. We collected the data for all human known 
protein coding genes (consisting in a total of 22.257 
entries with a stable Ensembl Gene Identifier (ENSG)). 
For each gene we then retained for further analysis only 
the longest Transcript Identifier (ENST). The data on 
human miRNAs were extracted from Ensembl v.57, that 
includes miRBase v. 13 (Table SI). To identify the iMSLs, 
we used eight tools for miRNA/target gene interaction 
predictions: TargetScan human v. 5.0 [3,82], miRanda - 
release 2008 [83,84], RNA22 [85], PITA-4way [86], Mir- 
Target2 [87], PicTar [88], Diana microT v.3 [89] and 
TargetMiner v.l [90]. To test the over- representation of 
the putative iMSLs, we performed two different types 
of randomization strategies. Specifically, we randomly 
permuted 1000 times the intronic miRNA/host-gene link 
and the union of miRNA/target gene interactions datasets 
predicted by the different algorithms. In both cases we 
created, according to the two reshuffling strategies, 1000 
independent reshuffled copies of the original network. 
Then for each of them we evaluated the number of iMSLs 
confirmed by at least one algorithm and obtained in this 
way the two histograms plotted in Figure 7C. Then for 
both reshuffling strategies a Z score can be defined as : 
Z = x ~< x> r t w here x is the actual number of iMSLs in 
the network (i.e. 77), while < x > r and o r are the mean 
value and the standard deviation of the distribution of 
iMSLs in the reshuffled samples. These Z scores turned 
out to be rather large: Z = 4.63 for the first test and 
Z = 5.52 for the second one (see Figure 7C). These values 



show that the number of iMSLs in the network is defi- 
nitely larger than random and thus suggest, in agreement 
with what already observed in [20,26], that this particular 
class of network motifs is most probably under positive 
selection. 

Additional file 



Additional file 1 : Supplementary information. Single pdf file containing 
the details of our modeling strategy, a comparison with alternative models 
of microRNA regulation, and the results of our bioinformatic analysis. 



Competing interests 

The authors declare that they have no competing interests. 
Author's contributions 

CB, MO and MC designed research. CB, MO, ME, DC performed research. CB, 
MO, ME and MC wrote the paper. All authors read and approved the final 
manuscript. 

Acknowledgements 

We would like to thank Marco Cosentino Lagomarsino for useful comments 
on the manuscript. MO acknowledges support from the Human Frontier 
Science Program Organization (Grant RGY0069/2009-C). DC acknowledges 
support by Fondazione Piemontese per la Ricerca sul Cancro-ONLUS and 
Associazione Italiana per la Ricerca sul Cancro (Intramural Grant 5x1 000 2008). 

Author details 

1 Human Genetics Foudation (HuGeF), V. Nizza 52, 1-10126, Torino, Italy. 2 Center 
for Complex Systems in Molecular Biology and Medicine, University of Torino, 
V. Accademia Albertina 1 3, 1-1 01 00 Torino, Italy. 3 Genomic Physics Group, UMR 
7238 CNRS "Microorganism Genomics", 75006, Paris, France. 4 Universite Pierre 
et Marie Curie, 1 5 rue de L'Ecole de Medecine, 75006, Paris, France. 5 National 
Research Council (CNR), Institute of Informatics and Telematics (I IT") and 
Institute of Clinical Physiology (IFC), Laboratory for Integrative System 
Medicine (LISM), Via Giuseppe Moruzzi 1, 1-56124, Pisa, Italy. 6 IRC@C: Institute 
for Cancer Research at Candiolo, School of Medicine, University of Torino, Str. 
Prov. 142, Km. 3.95, Candiolo 1-1 0060, Torino, Italy. 7 Dipartimento di Fisica 
Teorica and INFN, University of Torino, V. Pietro Giuria 1, 1-1 01 25, Torino, Italy. 

Received: 1 1 June 201 2 Accepted: 1 October 201 2 
Published: 1 0 October 201 2 

References 

1 . Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked 
by adenosines, indicates that thousands of human genes are 
microRNA targets. Cell 2005, 1 20:1 5-20. [http://dx.doi.Org/1 0.1 01 6/j.cell. 
2004.12.035] 

2. Flynt AS, Lai EC: Biological principles of microRNA-mediated 
regulation: shared themes amid diversity. Nat Rev Genet 2008, 
9(1 1):83 1-842. [http://dx.doi.org/10.1038/nrg2455] 

3. Friedman RC, Farh KKH, Burge CB, Bartel DP: Most mammalian mRNAs 
are conserved targets of microRNAs. Genome Res 2009, 1 9:92-1 05. 
[http://dx.doi.Org/1 0.1 101/g r.08270 1.108] 

4. Ambros V: The functions of animal microRNAs. Nature 2004, 
431(7006):350-355. [http://dx.doi.org/10.1038/nature02871] 

5. Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and 
function. Cell 2004, 1 1 6(2):28 1 -297. 

6. Bushati N, Cohen SM: microRNA functions. Annu Rev Cell Dev Biol 2007, 
23:175-205. [http://dx.doi.org/10.1146/annurev.cellbio.23.090506. 
123406] 

7. Stefani G, Slack FJ: Small non-coding RNAs in animal development. Nat 

Rev Mol Cell Biol 2008, 9(3):21 9-230. [http://dx.doi.org/10.1038/nrm2347] 

8. Alvarez-Garcia I, Miska EA: MicroRNA functions in animal development 
and human disease. Development 2005, 1 32(21 ):4653-4662. [http://dx. 
doi.org/10.1242/dev.02073] 



Bosia etal. BMC Systems Biology 201 2, 6:1 31 
http://www.biomedcentral.eom/1 752-0509/6/1 31 



9. Esquela-Kerscher A, Slack FJ: Oncomirs - microRNAs with a role in 

cancer. Nat Rev Cancer 2006, 6(4):259-269. [http://dx.doi.Org/1 0.1 038/ 
nrc1840] 

1 0. Hinske LCG, Galante PAF, Kuo WP, Ohno-Machado L: A potential role for 
intragenic miRNAs on their hosts' interactome. BMC Genomics 201 0, 

1 1 :533. [http://dx.doi.Org/1 0.1 1 86/1 471 -21 64- 1 1 -533] 

11. Li SC, Tang P, Lin WC: Intronic microRNA: discovery and biological 
implications. DNA Cell Biol 2007, 26(4):1 95-207. [http://dx.doi.Org/1 0. 
1089/dna.2006.0558] 

1 2. Saini HK, Griffiths-Jones S, Enright AJ: Genomic analysis of human 
microRNA transcripts. Proc Natl Acad Sci U S A 2007, 

1 04(45):1 771 9-1 7724. [http://dx.doi.Org/1 0.1 073/pnas.07038901 04] 

1 3. Rodriguez A, Griffiths-Jones S, Ashurst JL, Bradley A: Identification of 
mammalian microRNA host genes and transcription units. Genome 
Res 2004, 14(1 0A):1 902-1 910. [http://dx.doi.org/10.1101/gr.2722704] 

14. Baskerville S, Bartel DP: Microarray profiling of microRNAs reveals 
frequent coexpression with neighboring miRNAs and host genes. 
RNA 2005, 1 1 (3):241 -247. [http://dx.doi.Org/1 0.1 261 /rna.7240905] 

1 5. Lin SL, Miller JD, Ying SY: Intronic microRNA (miRNA). J Biomed 
Biotechnol 2006, 2006(4):2681 8. [http://dx.doi.Org/1 0.1 1 55/JBB/2006/ 
26818] 

1 6. Kim YK, Kim VN: Processing of intronic microRNAs. EMBO J 2007, 
26(3)775-783. [http://dx.doi.Org/1 0.1 038/sj.emboj.7601 51 2] 

1 7. Ma N, Wang X, Qiao Y, Li F, Hui Y, Zou C, Jin J, Lv G, Peng Y, Wang L, 
Huang H, Zhou L, Zheng X, Gao X: Coexpression of an intronic 
microRNA and its host gene reveals a potential role for miR-483-5p 
as an IGF2 partner. Mol Cell Endocrinol 201 1 , 333:96-1 01 . [http://dx.doi. 
org/1 0.1 01 6/j.mce.201 0.1 1.027] 

18. Kim VN, Han J, Siomi MC: Biogenesis of small RNAs in animals. Nat Rev 
Mol Cell Biol 2009, 1 0(2):1 26-1 39. [http://dx.doi.Org/1 0.1 038/nrm2632] 

1 9. Morlando M, Ballarino M, Gromak N, Pagano F, Bozzoni I, Proudfoot NJ: 
Primary microRNA transcripts are processed co-transcriptionally. 
Nat Struct Mol Biol 2008, 1 5(9):902-909. 

20. Ozsolak F, Poling LL, Wang Z, Liu H, Liu XS, Roeder RG, Zhang X, Song JS, 
Fisher DE: Chromatin structure analyses identify miRNA promoters. 
Genes Dev 2008, 22(22):31 72-31 83. [http://dx.doi.Org/1 0.1 1 01 /gad. 
1706508] 

21 . Monteys AM, Spengler RM, Wan J, Tecedor L, Lennox KA, Xing Y, Davidson 
BL: Structure and activity of putative intronic miRNA promoters. RNA 
201 0, 1 6(3):495-505. [http://dx.doi.Org/1 0.1 261 /rna.1 73 1 91 0] 

22. He C, Li Z, Chen P, Huang H, Hurst LD, Chen J: Young intragenic miRNAs 
are less coexpressed with host genes than old ones: implications of 
miRNA-host gene coevolution. Nucleic Acids Res 201 2, 40(9):4002-401 2. 
[http://dx.doi.Org/1 0.1 093/nar/gkr1 31 2] 

23. Barik S: An intronic microRNA silences genes that are functionally 
antagonistic to its host gene. Nucleic Acids Res 2008, 36(1 6):5232-5241 . 
[http://dx.doi.org/10.1093/nar/gkn513] 

24. Lutter D, Marr C, Krumsiek J, Lang EW: Intronic microRNAs support 
their host genes by mediating synergistic and antagonistic 
regulatory effects. BMC Genomics 2010, 1 1:224. [http://dx.doi.org/10. 
1186/1471-2164-11-224] 

25. Tsang J, Zhu J, van Oudenaarden A: MicroRNA-mediated feedback and 
feedforward loops are recurrent network motifs in mammals. Mol 
Cell 2007, 26(5)753-767. [http://dx.doi.Org/1 0.1 01 6/j.molcel.2007.05.01 8] 

26. Megraw M, Sethupathy P, Gumireddy K, Jensen S, Huang Q: Isoform 
specific gene auto-regulation via miRNAs: a case study on miR-1 28b 
and ARPP-21 . Theor Chem Acc 2009, 1 25:593-598. 

27. Alon U: Network motifs: theory and experimental approaches. Nat 
Rev Genet 2007, 8(6):450-461 . [http://dx.doi.Org/1 0.1 038/nrg21 02] 

28. Sun Y, Bai Y, Zhang F, Wang Y, Guo Y: miR-1 26 inhibits non-small cell 
lung cancer cells proliferation by targeting EGFL7. Biochem Biophys 
Res Commun 201 0, 391 (3):1 483-1 489. [http://dx.doi.Org/1 0.1 01 6/j.bbrc. 
2009.12.098] 

29. Nikolic I, Plate KH, Schmidt MH: EGFL7 meets miRNA- 1 26: an 
angiogenesis alliance. J Angiogenes Res 201 0, 2:9. [http://dx.doi.Org/1 0. 
1186/2040-2384-2-9] 

30. Alon U: An Introduction to Systems Biology: Design Principles of Biological 
Circuits. Boca Raton; FL: Chapman & Hall/CRC Mathematical^aylor &amp 
Francis; 2007. 



Page 15 of 16 



31 . Inui M, Martello G, Piccolo S: MicroRNA control of signal transduction. 

Nat Rev Mol Cell Biol 201 0, 1 1 (4):252-263. [http://dx.doi.Org/1 0.1 038/ 
nrm2868] 

32. Li X, Cassidy JJ, Reinke CA, Fischboeck S, Carthew RW: A microRNA 
imparts robustness against environmental fluctuation during 
development. Cell 2009, 1 37(2):273-282. [http://dx.doi.Org/1 0.1 01 6/j. 
cell.2009.01.058] 

33. Li Y, Wang F, Lee JA: MicroRNA-9a ensures the precise specification of 
sensory organ precursors in Drosophila. Genes Dev 2006, 
2O(20):2793-2805. [http://dx.doi.Org/10.1 101 /gad.1 466306] 

34. Ebert MS, Sharp PA: Roles for MicroRNAs in Conferring Robustness to 
Biological Processes. Cell 2012, 149(3):5 15-524. [http://dx.doi.org/10. 

1 01 6/j.cell.201 2.04.005] 

35. Haynes KA, Silver PA: Eukaryotic systems broaden the scope of 
synthetic biology. J Cell Biol 2009, 187(5):589-596. [http://dx.doi.org/10. 
1083/jcb.200908138] 

36. Khalil AS, Collins JJ: Synthetic biology: applications come of age. Nat 
Rev Genet 201 0, 1 1 (5):367-379. [http://dx.doi.Org/1 0.1 038/nrg2775] 

37. Deans TL, Cantor CR, Collins JJ: A tunable genetic switch based on 
RNAi and repressor proteins for regulating gene expression in 
mammalian cells. Cell 2007, 1 30(2):363-372. [http://dx.doi.Org/1 0.1 01 6/ 
j.cell.2007.05.045] 

38. Rinaudo K, Bleris L, Maddamsetti R, Subramanian S, Weiss R, Benenson Y: 
A universal RNAi-based logic evaluator that operates in mammalian 
cells. Nat Biotechnol 2007, 25(7)795-801 . [http://dx.doi.Org/1 0.1 038/ 
nbt1307] 

39. Bleris L, Xie Z, Glass D, Adadey A, Sontag E, Benenson Y: Synthetic 
incoherent feedforward circuits show adaptation to the amount of 
their genetic template. Mol Syst Biol 201 1 , 7:51 9. [http://dx.doi.Org/1 0. 
1038/msb.2011.49] 

40. Bintu L, Buchler NE, Garcia HG, Gerland U, Hwa T, Kondev J, Phillips R: 
Transcriptional regulation by the numbers: models. CurrOpin Genet 
Dev 2005, 1 5(2):1 1 6-1 24. [http://dx.doi.Org/1 0.1 01 6/j.gde.2005.02.007] 

41. Valencia-Sanchez MA, Liu J, Hannon GJ, Parker R: Control of translation 
and mRNA degradation by miRNAs and siRNAs. Genes Dev 2006, 
20(5):5 15-524. [http://dx.doi.Org/10.1 101 /gad.1 399806] 

42. Cuccato G, Polynikis A, Siciliano V, Graziano M, di Bernardo M, di Bernardo 
D: Modeling RNA interference in mammalian cells. BMC Syst Biol 201 1, 
5:1 9. [http://dx.doi.Org/1 0.1 1 86/1 752-0509-5- 1 9] 

43. Komorowski M, Miekisz J, Kierzek AM: Translational repression 
contributes greater noise to gene expression than transcriptional 
repression. Biophys J 2009, 96(2):372-384. [http://dx.doi.Org/10.1016/j. 
bpj.2008.09.052] 

44. Osella M, Bosia C, Cora D, Caselle M: The role of incoherent microRNA- 
mediated feedforward loops in noise buffering. PLoS Comput Biol 
201 1 , 7(3):e1 001 1 01 . [http://dx.doi.Org/1 0.1 371 /journal.pcbi.1 001 101] 

45. Levine E, Zhang Z, Kuhlman T, Hwa T: Quantitative characteristics of 
gene regulation by small RNA. PLoSBiol 2007, 5(9):e229. [http://dx.doi. 
org/1 0.1 371/journal.pbio.0050229] 

46. Shimoni Y, Friedlander G, Hetzroni G, Niv G, Altuvia S, Biham O, Margalit H: 
Regulation of gene expression by small non-coding RNAs: a 
quantitative view. Mol Syst Biol 2007, 3:1 38. [http://dx.doi.Org/1 0.1 038/ 
msb4100181] 

47. Mehta P, Goyal S, Wingreen NS: A quantitative comparison of 
sRNA-based and protein-based gene regulation. Mol Syst Biol 2008, 
4:221. [http://dx.doi.org/10.1038/msb.2008.58] 

48. Mukherji S, Ebert MS, Zheng GXY, Tsang JS, Sharp PA, van Oudenaarden 
A: MicroRNAs can generate thresholds in target gene expression. 
Nat Genet 201 1, 43(9):854-859. [http://dx.doi.org/10.1038/ng.905] 

49. Krol J, Loedige I, Filipowicz W: The widespread regulation of microRNA 
biogenesis, function and decay. Nat Rev Genet 201 0, 1 1 (9):597-61 0. 
[http://dx.doi.org/10.1038/nrg2843] 

50. Tigges M, Denervaud N, Greber D, Stelling J, Fussenegger M: A synthetic 
low-frequency mammalian oscillator. Nucleic Acids Res 201 0, 
38(8):2702-271 1. [http://dx.doi.org/10.1093/nar/gkq121] 

51 . Kaern M, Elston TC, Blake WJ, Collins JJ: Stochasticity in gene 
expression: from theories to phenotypes. Nat Rev Genet 2005, 
6(6):45 1 -464. [http://dx.doi.Org/1 0.1 038/nrg 1615] 

52. Maheshri N, O'Shea EK: Living with noisy genes: how cells function 
reliably with inherent variability in gene expression. Annu Rev 



Bosia etal. BMC Systems Biology 201 2, 6:1 31 
http://www.biomedcentral.eom/1 752-0509/6/1 31 



Page 16 of 16 



Biophys Biomol Struct 2007, 36:413-434. [http://dx.doi.Org/10.1 146/ 
annurev.biophys.36.040306.1 32705] 

53. Raj A, van Oudenaarden A. Ce//2008, 135(2):21 6-226. [http://dx.doi.org/ 
1 0.1 01 6/j.cell.2008.09.050] 

54. Re A, Cora D, Taverna D, Caselle M: Genome-wide survey of microRNA- 
transcription factor feed-forward regulatory circuits in human. Mol 
Biosyst 2009, 5(8):854-867. [http://dx.doi.org/10.1039/b900177h] 

55. El Baroudi, M, Cora D, Bosia C, Osella M, Caselle M: A curated database of 
miRNA mediated feed-forward loops involving MYC as master 
regulator. PLoSOne 201 1, 6(3):e1 4742. [http://dx.doi.org/10.1371/ 
journal. pone.0014742] 

56. Eduati F, Camillo BD, Karbiener M, Scheideler M, Cora D, Caselle M, Toffolo 
G: Dynamic modeling of miRNA-mediated feed-forward loops. 

J Comput Biol 201 2, 1 9(2):1 88-1 99. [http://dx.doi.Org/1 0.1 089/cmb.201 1 . 
0274] 

57. Rosenfeld N, Elowitz MB, Alon U: Negative autoregulation speeds the 
response times of transcription networks. J Mol Biol 2002, 

323(5)785-793. 

58. Mangan S, Itzkovitz S, Zaslaver A, Alon U: The incoherent feed-forward 
loop accelerates the response-time of the gal system of Escherichia 
coli. J Mol Biol 2006, 356(5):1 073-1 081 . [http://dx.doi.Org/1 0.1 01 6/j.jmb. 
2005.12.003] 

59. Mangan S, Zaslaver A, Alon U: The coherent feedforward loop serves 
as a sign-sensitive delay element in transcription networks. J Mol Biol 
2003, 334(2):1 97-204. 

60. Fan J, Yang X, Wang W, Wood WH, Becker KG, Gorospe M: Global 
analysis of stress-regulated mRNA turnover by using cDNA arrays. 
Proc Natl Acad Sci USA 2002, 99(1 6):1 061 1 -1 061 6. [http://dx.doi.Org/1 0. 
1073/pnas.162212399] 

61 . Chen CYA, Ezzeddine N, Shyu AB: Messenger RNA half-life 
measurements in mammalian cells. Methods Enzymol 2008, 
448:335-357. [http://dx.doi.Org/1 0.1 01 6/S0076-6879(08)0261 7-7] 

62. Yang E, van Nimwegen, E, Zavolan M, Rajewsky N, Schroeder M, 
Magnasco M, Darnell JE: Decay rates of human mRNAs: correlation 
with functional characteristics and sequence attributes. Genome Res 
2003, 1 3(8):1 863-1 872. [http://dx.doi.Org/1 0.1 1 01 /gr.1 272403] 

63. Sharova LV, Sharov AA, Nedorezov T, Piao Y, Shaik N, Ko MSH: Database 
for mRNA half-life of 19 977 genes obtained by DNA microarray 
analysis of pluripotent and differentiating mouse embryonic stem 
cells. DNA Res 2009, 16:45-58. [http://dx.doi.org/10.1093/dnares/dsn030] 

64. Yen HCS, Xu Q, Chou DM, Zhao Z, Elledge SJ: Global protein stability 
profiling in mammalian cells. Science 2008, 322(5903):91 8-923. [http:// 
dx.doi.org/10.1 1 26/science.1 160489] 

65. van Rooij, E, Sutherland LB, Qi X, Richardson JA, Hill J, Olson EN: Control 
of stress-dependent cardiac growth and gene expression by a 
microRNA. Science 2007, 31 6(5824):575-579. [http://dx.doi.Org/1 0.1 1 26/ 
science.1 139089] 

66. Kai ZS, Pasquinelli AE: MicroRNA assassins: factors that regulate the 

disappearance of miRNAs. Nat Struct Mol Biol 20 1 0, 1 7:5- 1 0. [http://dx. 
doi.org/10.1038/nsmb.1762] 

67. Sethi P, Lukiw WJ: Micro-RNA abundance and stability in human 
brain: specific alterations in Alzheimer's disease temporal lobe 
neocortex. Neurosci Lett 2009, 459(2): 1 00-1 04. [http://dx.doi.Org/1 0. 
1016/j.neulet.2009.04.052] 

68. Chatterjee S, Grosshans H: Active turnover modulates mature 
microRNA activity in Caenorhabditis elegans. Nature 2009, 
461 (7263):546-549. [http://dx.doi.Org/1 0.1 038/nature08349] 

69. Kollmann M, Lovdok L, Bartholomew K, Timmer J, Sourjik V: Design 
principles of a bacterial signalling network. Nature 2005, 
438(7067):504-507. [http://dx.doi.Org/1 0.1 038/nature04228] 

70. Matthews HR, Reisert J: Calcium, the two-faced messenger of olfactory 
transduction and adaptation. CurrOpin Neurobiol 2003, 1 3(4):469-475. 

71 . Ma W, Trusina A, El-Samad H, Lim WA, Tang C: Defining network 
topologies that can achieve biochemical adaptation. Cell 2009, 
1 38(4)760-773. [http://dx.doi.Org/1 0.1 01 6/j.cell.2009.06.01 3] 

72. Sontag ED: Remarks on feedforward circuits, adaptation, and pulse 
memory. lETSyst Biol 201 0, 4:39-51 . [http://dx.doi.Org/1 0.1 049/iet-syb. 
2008.0171] 



73. Francois P, Siggia ED: A case study of evolutionary computation of 

biochemical adaptation. Phys Biol 2008, 5(2):026009. [http://dx.doi.org/ 
1 0.1 088/1 478-3975/5/2/026009] 

74. Goentoro L, Kirschner MW: Evidence that fold-change, and not 
absolute level, of beta-catenin dictates Wnt signaling. Mol Cell 2009, 
36(5):872-884. [http://dx.doi.Org/10.1016/j.molcel.2009.1 1 .01 7] 

75. Cohen-Saidon C, Cohen AA, Sigal A, Liron Y, Alon U: Dynamics and 
variability of ERK2 response to EGF in individual living cells. Mol Cell 
2009, 36(5):885-893. [http://dx.doi.Org/10.1016/j.molcel.2009.1 1.025] 

76. Goentoro L, Shoval O, Kirschner MW, Alon U: The incoherent 
feedforward loop can provide fold-change detection in gene 
regulation. Mol Cell 2009, 36(5):894-899. [http://dx.doi.Org/1 0. 1 0 1 6/j. 
molcel.2009.1 1.018] 

77. Dublanche Y, Michalodimitralis K, Kummerer N, Foglierini M, Serrano L: 
Noise in transcription negative feedback loops: simulation and 
experimental analysis. Molecular Systems Biology 2006, 2:41 . 

78. Shahrezaei V, Ollivier J, PS S: Colored extrinsic fluctuations and 
stochastic gene expression. Mol Syst Biol 2008, 4:1 96. 

79. Singh A, Hespanha JP: Optimal feedback strength for noise 
suppression in autoregulatory gene networks. Biophys J 2009, 
96(1 0):401 3-4023. [http://dx.doi.Org/1 0.1 01 6/j.bpj.2009.02.064] 

80. Baek D, Villen J, Shin C, Camargo FD, Gygi SP, Bartel DP: The impact of 
microRNAs on protein output. Nature 2008, 455(7209):64-71 . [http:// 
dx.doi.org/1 0.1 038/nature07242] 

81 . Arvey A, Larsson E, Sander C, Leslie CS, Marks DS: Target mRNA 
abundance dilutes microRNA and siRNA activity. Mol Syst Biol 201 0, 
6:363. [http://dx.doi.Org/1 0. 1 038/msb.20 1 0.24] 

82. Lewis BP, hung Shih, I, Jones-Rhoades MW, Bartel DP, Burge CB: 
Prediction of mammalian microRNA targets. Cell 2003, 
115(7)787-798. 

83. John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks: Human MicroRNA 
targets. PLoSBiol 2004, 2(1 1):e363. [http://dx.doi.org/10.1371/journal. 
pbio.0020363] 

84. Betel D, Wilson M, Gabow A, Marks DS, Sander C: The microRNA.org 
resource: targets and expression. Nucleic Acids Res 2008, 36(Database 
issue):D1 49-D1 53. [http://dx.doi.Org/1 0.1 093/nar/gkm995] 

85. Miranda KC, Huynh T, Tay Y, Ang YS, Tarn WL, Thomson AM, Lim B, 
Rigoutsos I: A pattern-based method for the identification of 
MicroRNA binding sites and their corresponding heteroduplexes. 
Cell 2006, 1 26(6):1 203-1 21 7. [http://dx.doi.Org/1 0.1 01 6/j.cell.2006.07.031] 

86. Kertesz M, lovino N, Unnerstall U, Gaul U, Segal E:The role of site 
accessibility in microRNA target recognition. Nat Genet 2007, 
39(1 0):1 278-1 284. [http://dx.doi.Org/1 0.1 038/ng21 35] 

87. Wang X, Naqa IME: Prediction of both conserved and nonconserved 
microRNA targets in animals. Bioinformatics 2008, 24(3):325-332. 
[http://dx.doi.Org/1 0.1 093/bioinformatics/btm595] 

88. Krek A, Grun D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, 
da Piedade I, Gunsalus KC, Stoffel M, Rajewsky N: Combinatorial 
microRNA target predictions. Nat Genet 2005, 37(5):495-500. [http://dx. 
doi.org/10.1038/ng1536] 

89. Maragkakis M, Reczko M, Simossis VA, Alexiou P, Papadopoulos GL, 
Dalamagas T, Giannopoulos G, Goumas G, Koukis E, Kourtis K, Vergoulis T, 
Koziris N, Sellis T, Tsanakas P, Hatzigeorgiou AG: DIANA-microT web 
server: elucidating microRNA functions through target prediction. 
Nucleic Acids Res 2009, 37(Web Server issue):W273— W276. [http://dx.doi. 
org/1 0.1 093/nar/gkp292] 

90. Bandyopadhyay S, Mitra R: TargetMiner: microRNA target prediction 
with systematic identification of tissue-specific negative examples. 
Bioinformatics 2009, 25(20):2625-263 1 . [http://dx.doi.Org/1 0. 1 093/ 
bioinformatics/btp503] 

91. Dekel E, Alon U: Optimality and evolutionary tuning of the 
expression level of a protein. Nature 2005, 436(7050):588-592. [http:// 
dx.doi.org/1 0.1 038/nature03842] 

92. Gillespie D: A general method for numerically simulating the 
stochastic time evolution of coupled chemical reactions. J Comput 
Phys 1976,22:403-434. 

/ \ 
doi:1 0.1 1 86/1 752-0509-6-1 31 

Cite this article as: Bosia etal.: Gene autoregulation via intronic microRNAs 
and its functions. BMC Systems Biology 201 2 6:1 31 . 
\ J 



